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It is shown that the theory of non-self -ad joint linear ordinary differential equations can 
be simplified and unified if, instead of specifying linear boundary conditions in the conven- 
tional way, one merely specifies the linear subspace determined by the boundary conditions. 
If this is done, the corresponding linear space of the adjoint problem is the orthogonal 
complement of the original space with respect to the scalar product defined by the right-hand 
side of Green's bound arv formula. 



1. Introduction 

The theory of boundary- value problems for or- 
dinary differential equations of the form 



*=0 



=/(*) 



(1) 



is now regarded as more or less complete. The self- 
adjoint case is presented in most textbooks owing 
to its importance in applications. A complete 
treatment of the non-self-adjoint problem is less 
common but may be found in Coddington, Levinson, 2 
and Ince. 3 However, in these as well as in the 
original publications, e.g., Bocher, 4 the boundary 
conditions are stated in the form 

ov«a), . . ., u (n - 1} (a),u(b), . • ., u {n - l) (&)) 

G*=l, ...,m) (2) 

where the co M are m linearly independent linear 
functions in the 2n dimensional number space 
R 2n (0^m^2ri). 

The m linear equations o> M =0 (m=1j • • •> m ) 
define a (2n— m) dimensional linear subspace 12 £ R 2n . 
The set of solutions u(t) of the problem (l)/(2) 
remains unchanged if the functions co M are replaced 
by another system of equations co M =0 (/a=1, . . .,ra) 
which define the same subspace 12 £ R 2n . 



1 The work was carried out while both authors were employed at the University 
of Maryland , the first in the Department of Mathematics, the second in the Insti- 
tute for Fluid Dynamics and Applied Mathematics. The paper was prepared 
while both authors were affiliated with the National Bureau of Standards, the 
first working under a National Bureau of Standards contract with The American 
University. The authors' present addresses are, respectively: Mathematics 
Institute, University of Zurich, Zurich, Switzerland; and Computing Center, 
Syracuse I'uiversify, Syracuse, N.Y. 

2 Earl A. Coddhmton, Norman Levinson, Theory of ordinary differential 
equations (McGraw-Hill Book Co., Inc., New York, N.Y., 1955). 

3 E. L. Ince, Ordinary differential equations, London (1927). 

4 M. Bocher, Amplications and generalizations of the conception of adjoint 
systems, Trans. Am. Math. Soc. 14, 403 to 420 (1913). 



This fact causes an unnecessary complication 
for the theoretical investigation of such boundary 
value problems. It results especially in a somewhat 
artificial definition of the adjoint problem. These 
difficulties do not occur if one prescribes instead of 
the boundary conditions (2) only the subspace 
12 of R 2n , disregarding the manner in which 12 is 
represented. If this is done, the corresponding 
linear space of the adjoint problem is the orthogonal 
complement of the original space 12 with respect to a 
scalar product defined by the right-hand side of 
Green's boundary value formula. 

It is shown in this article that the theory of the 
non-self-adjoint problems can be simplified and 
unified by a consequent use of these ideas and the 
methods of linear algebra. 

2. Green's Formula 

As usual we denote by C* (k=Q,l, . . .) the space 
of all real-valued 5 funetions u(t) which possess k 
continuous derivatives on a given closed interval 
[a, b]. Let L[u] be the linear differential operator of 
order n(^l) defined by 



L[u] 



*2>.-»(0« (, "«) 

v=l 



(3) 



where a n - v (t) £ C v ,u(t) £ C n and a (t)^0 for [a, b]. 
To every function u(t) £ C n ~ l one can compute the 
2n real numbers 



u(a) ,u' (a) , . . ., u in - 1] (a) ,u(b) , u' (b) , . . 



#(«-!) 



(b). (4) 



We regard these 2n numbers as the components of a 
vector u of the 2n-dimensional linear number space 
R 2n and call it the boundary value vector of u(t) . 
The correspondence between the function u(t) 
and its boundary value vector u is a linear mapping 

5 The restriction to real-valued functions is not essential, as can be readily 
ascertained. 
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u= P (u) 



(u(t)£ C n ~\ u£R 2n ) 



(5) 



from C n ~ l onto R 2n . 

Let 12 be a given subspace of the R 2n . Then the 
boundary value problem in consideration assumes 
the form 



L[u]=f(t), p(u)£Q, 



(*n) 



with given j(t) £ C°. The corresponding homoge- 
neous problem is 

L[u] = 0, p(u)£Q. (t h ) 

For any two functions u(t), v(t) £ C n , the integral 

(L[u],v)= f b L[u]vdt (6) 

J a 

is defined and by repeated partial integrations can 
be written as 



where 



(L[u], v) = (u, L*[v]) + [u,v]l, 



£n*]=S(-l)'(a*-,) w , 

r=0 



M=X S (-l)"-'-H* w (o i ,- F t>)''-'- 1 |». 

y = /X=f+1 



(7) 
(8) 
(9) 



Equation (7) is the well-known Green's formula. 

The linear differential operator L* is called the 
adjoint operator of L. The expression [u,v]l is a 
bilinear form in the two boundary value vectors 
p(u), p(v). For arbitrary vectors 

x={.., £«>,..;. .,£»>,..}, y={.., ,<«>,. .;..,!!»),..} 

of the R 2n , this bilinear form can be written in the 
form 



P 



W — 1 ft — v— 1 ?1 — J>— 1 /\\ 

= 2 S S(-D X Q) 



y=0 m=0 

•K_-A-i(&)^ , ^ 4 '-^ x --A-i(«)^ a> '7i < "] 

With this notation we get Green's formula, 
(Z[«L0) = (i*,£*M)+*(p(«),p(t>)). 



(10) 



(11) 



It is obvious that the correspondence L — > L* is 
linear. Furthermore, the adjoint operator to L* is 
again L. To show this, apply Green's formula to 
any u,v £ C re with p(u)=p(v) =0: 

(L[u],v) -(u,L*M) =0 

(Z*M, «)-(», £**M)=0, 

which results in 

(XM-z**m,v)=.o 



(for all w,#(: C w with p(u) = p(v) = 0) and thus leads 
to 

L[u]=L**[u] 

for all functions u(t) £ C n . In the last step one has 
to use the same technique used in calculus of varia- 
tions for the proof of Euler's condition. 

The differential operator L is said to be self -ad joint 
if Z*=Z. On comparing the forms (3) and (8) of 
L and Z* one sees immediately that n must be even 
if the operator is self-adjoint. 

In the case of a self -ad joint operator L the left- 
hand side of (11) is skew-symmetric with respect to 
u and v, and this implies that the bilinear form 
<p(u,\) is also skew-symmetric. 

We will now show that <£>(x,y) is not degenerate. 
To prove this it is sufficient to verify the following 
statement: If for a fixed vector y£ R 2n we have 
<p(x,y) = for every x£ R 2n , then it follows that 

y=o. 

We prove this first for a fixed vector y of the form 



y«=M 8> , 






. Vn-U 



..,0} 



or 



y 6 ={0, . . ., 0, r\ 



(6) 



'3 Vn-l 



(&) A 



(12a) 



(12b) 



taking x correspondingly as x a or x 6 , analogously de- 
fined. In this case the expression <p(x a , y a ) for a=a 
or a = b is a bilinear form in the two ^-dimensional 
vectors {^\ . . ., &l\} and {Vo a) , . . ., V A}. 
Therefore, if <p(x a , J a )=0 holds for every x a and a 
fixed j^ then in expression (10) the coefficient of 
every £ ( J° must vanish separately. Thus 

*S' "I? (-l) x ft) aftrA-i («) v^=0 

(v=0, . . .,n-l). 

For v=n—l this has the form a (a) v ( o ) = } and since 
a (a) 9^0 it follows that 



vV-- 



=0. 



For ?7—7i—2 we get now 
and 



-a (a) rj^ = 



Continuing in this way, we ultimately show that 

vf=0 (m=0, . . ., n-l) 

and thus y«=0. 

Now the proof is readily completed. From defini- 
tion (10) it is seen that for arbitrary x=x a +x 6 and 

y=y«+y &J 

<p(x, jr)=«p(x 6 , y 6 )— <p(x a , Jo). 



84 



For a fixed y this expression si 1 all be zero for every 
x£R 2n and in particular for arbitrary vectors x=x fl 
or x=x 6 . Hence the two terms <p(x a , y a ) and <p(x b , j b ) 
have to vanish separately for every x a and x 6 , vr- 
spectively. In accordance with the first part of the 
proof, this is only possible for y a =0 and y&=0, i.e., 
y=0. 

3. The Adjoint Problem 

The bilinear form <p was shown to be nondegener- 
ate. Accordingly, it can be employed to define 
"orthogonality" of vectors x, y£R 2n . We will say 
x and y are ^-orthogonal if <p(x, y)=0. It should, be 
observed that this ^-orthogonality is not necessarily 
symmetrical. 

One can speak, in particular, of the ^-orthogonal 
space to 12 that is the space 1 - of all vectors y£R 2n 
which are ^-orthogonal to all vectors x£l2. 

For the dimension of 12- 1 evidently holds 



dim 12- L = 2ft— dim 1] = 



-m. 



(13) 



Equation (13) implies 

dim (Q x ) J -=dim U, 
and hence, together with (ft 1 )- 1 ^, that 

(oj-)J-=o. 



(14) 



It is worth while pointing out that the relation 
^n^" 1 "^^ is not, in general, valid; there may exist a 
nonzero vector x£l2 which is ^-orthogonal to all 
vectors of U. 

With this terminology we define the adjoint 
boundary value problem of (ir N ) : 

(n*) L*[u]=g®, p(u)£Q 

with given g(t) ^ 0°. 

From L**=L and (14), it follows that the adjoint 
problem of (t%) is again (tt n ) . 

In particular, (w N ) will be termed self-adjoint 6 if 



4. Homogeneous Problems 



(15) 



Consider the homogeneous problems (w H ) and 
(irff). It is known from the theory of homogeneous 
linear differential equations that all solutions u ^C n 
of the equation 

L[u]=0 (16) 

form a n-dimensional function space K. The solu- 
tions v£C n of L*[v] = form another n -dimensional 
function space K*. These two spaces K and K* 
are mapped by the operator p into two subspaces 

K= P (K) and K* = p(K*) of R 2 \ 



6 For self-adjoint problems Green's formula states 

(L[u], v) = (u,L[v]) for p(u), P (v) e O 

Usually one uses this relation to define self-adjointness of (ttn). The two defini- 
tions are equivalent as can be shown in a way similar to that used for the proof of 
L*=L. 



In these subspaces K and K* the operator p lias a 
unique inverse p -1 ; hence, K=p~ l (K) and K*= 
p -1 (K*). This is immediately clear because p(u) = () 
implies that all initial values u {v) (a) (p=0,1> . . ., 
n— 1) are zero; hence u(t) = owing to the unique- 
ness theorem for the initial value problem of eq (16). 

The operator L and the bilinear form <p(x,y) deter- 
mine the ^-orthogonal space K = p(K) A - of K. 

Theorem 1: The subspaces p(K*) and K A - = p(K)- L 
of R 2n are identical, i.e., 



K-l= p (K*). 



(17) 



Pkoof: For any two functions u(t) £K, v(t) £K* 
Green's formula gives <p(p(u), p(i>)) = 0, which implies 



p(K*)£ P (K)J-=K\ 



(18) 



Since the mapping p is regular in K and K*, it follows 
on the other hand that 



and 



Thus, we find 



dim p(jK)=dim K=n 
dim p(2f*)=dim K*=n. 

dim K ± =2n— dim K=n, 



and the inclusion (18) must constitute equality. 

To this point the boundary conditions [in (tt h ) and 
(t%) have not been used. To consider them now we 
introduce the space A of all admissible functions of 
(t h ), i.e., the space of all functions u(t)^C n with 
p(u)eti. Thus, the solutions of (ir H ) are, exactly the 
functions 

u(t)£A[\K. (19a) 

For these functions the boundary value vectors are 

p(u)^Sl(\p{K). (19b) 

We handle (ir%) in the same way. Let A* be the 
space of all v(t) £ C n for which p(u) £ 12. Then, the 
solutions of (tth) are exactly the functions 



v(f)£A*riK*. 



(20a) 



In view of relation (17) the boundary value vectors 
of these functions are 



P (v)^Q- L r\K^=^n P (K)\ 



(20b) 



Thus, the number of linearly independent solutions 
of (w H ) or (th) is given by 



i=dim (A()K)=dim (Sin P (K)) 



(21) 



and 



&*=dim (A()K*)=dim (^DpOfiQ-L), (22) 
respectively. 
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Theorem 2: The numbers k and k* of linearly 
independent solutions of (ir H ) and (r£) are connected 
by the relation 



] c * = m —n J rk. 



(23) 



To prove this we use the following well-known 
lemma from linear algebra: 

Let A be a linear space and A u A 2 two subspaces, 
then 

dim Ai+dim ^ 2 =dim(An^4 2 )+dim (A 1 +A 2 ), 

(24) 
where Ai-}-A 2 is the linear closure of A x and A 2 . 

This relation is employed in a slightly different 
form. We introduce the orthogonal subspaces 
Aj-^A and A\£_A of A l and A 2 with respect to 
a given (nondegenerate) bilinear form. Whence, 
we have 

dim A\=n— dim A± 
and 

dim (4irU£)=dim (A 1 +A 2 ) ± =n- dim (A^A 2 ). 

Thus (24) becomes 

dim A 2 — dim (A^A^ =dim A\ — dim (A\[\Ai). 

(25) 
In the present case, 

A=R 2n , A x = ^, A 2 = P (K); 

therefore 

dim A 1 =2n—m, dim ^4i = m, dim A 2 =n; 
hence, 

n-dim (nf]p(K))=m- dim (^nplif) 1 ), 
which completes the proof. 

5. The Nonhomogeneous Problem 

Consider the nonhomogeneous problem (t n ) to- 
gether with its adjoint homogeneous problem (th). 

We shall prove the ''alternative" theorem of the 
theory of systems of linear equations. 

Theorem 3 (Alternative theorem): (t n ) has a solu- 
tion if and only if (f,v)=0 for all solutions v(t) of 
(tS) . This can be expressed in the form 7 



L[A] = (A*(]K*)K 



(26) 



Here the symbol Jj_ denotes orthogonality in the 
function space C° with respect to the scalar product 

if, 9)- 

i L[A] is the subspace of C° containing all functions L[u] with u£A. 



Proof: (1) Assume (t&) is solvable and u(t) is a 
solution. Then, for any solution v(t) of (t|), it fol- 
lows from Green's formula that 

(f,v) = (L[u],v) = [u,L*[v]) +<p{ P {u) , P (»)). 

Since L*[v] = 0, then (u,L*[v])=0; and from p{u) £ £2 
and p(v) (^Q 1 - it follows that <p(p(u), p(v))=Q. Thus 
it is seen that (f,v) = 0. 

(2) Assume that (f,v) = for all solutions of (th). 
We shall show that (w N ) has a solution. 

From the existence theorem for ordinary differential 
equations it is known that the differential equation 
L[u]=f must have at least one solution u (t). 
For this function u (t) and any solution v(t) of (wh) 
we find from Green's formula that 

= (f, v) = (L[uo], ») = («•, L*[v])+<p(p(iio), p(t») 

=<p(p(u ), P (v)). 

Thus p(u ) is ^-orthogonal to all p(v) ^^CipiK) 1 - 
and 

p(u Q )e(^np(K)^) = tt + p(K) 

where the + sign again denotes the linear closure 
of 12 and p(K) . 

This inclusion states the existence of two vectors 
ill £ &, u 2 £ p(K) such that p(u ) = U!+u 2 . From the 
fact that p has a unique inverse p _1 in K it follows 
that there is exactly one function u 2 {t) = p~ l (vl 2 ) G^ 
for which L[u 2 } = and p(u 2 )=u 2 . We form a new 
function u(f)=u (t)—u 2 (f). This function u(t) is a 
solution of (ir N ), for obviously L[u]=f and p(u) = 
p(u ) — p(u 2 ) = n 1 ^Q. 

Taking the orthogonal complement on both sides 
of (26) we find 



l[a\ jl = (A*nK*)^=A*r\K*. 

Thus we get the following 

Corollary: If for a given function v(t) £ C°, 



(27) 



(v(t),L[u])=Q for every u(t)£A, 
then L*[v] = Q, p(v) £ Q and in particular v(t) £ C n . 

6. Green's Function 

The solution of the nonhomogeneous problem 
(ttn) — provided that it is solvable — is determined 
up to a solution of the homogeneous problem (ir H ). 
In order to obtain an explicit formula for the solution 
of (ir N ), it is necessary to assume an additional 
condition which uniquely determines the solution 
of (ir N ). Furthermore, it is desirable to modify (ir N ) 
in such a way that it is always solvable. 

With these two aims in mind we introduce an 
orthonormalized basis in each of the spaces 

A[\K and A*()K*. 
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Let 

u u u 2 , . . ., u k and v lf v 2 , . . ., v k (28) 
be these two bases. Witli the functions v v (t) we 
define the new boundary value problem 

£N=/(0-S (f(t)Mt))v v (t), p{u) e o. (29) 

Owing to the orthonormality of the v v (t) , the right- 
hand side of the differential equation is orthogonal 
to all v v (t) (v=l, . . ., k), and thus to all solu- 
tions of (V*/). Therefore (29) is solvable for every 
function f(t) £ C°. To find an additional condition 
which uniquely determines the solution of problem 
(29) , we recall the following fact: If u(t) is a solution 
of (29), then 



is also a solution for any real coefficients y m . 
can determine the y m in such a way that 



We 






{v=\, . . .,k). 



Owing to the orthonormality of the u^ these rela- 
tions are satisfied for 



Tm= 



-—(u(t)Mt)) (m=i, 



.,*). 



This shows that there is always a solution of (29), 
which is orthogonal to all u^ and thus to the whole 
space AV\K; obviously, this solution is uniquely 
determined. 

Consider instead of (t n ) the following problem: 



L[«]=/(0-S U,%) %(t) 

K=l 



p(*)€o 



fav) 



o 



(M=l, 



.,*). 



We know that this problem (£#) has exactly one 
solution u(t) £ C n for every function f(t) £ 0°. The 
correspondence between u{t) and j(t) is denoted by 
an operator G, i.e., 



u(t) = G[f(t)}. 



(30) 



This operator is obviously linear. It will be shown 
that G is an integral operator of the form 



/: 



Mt) = OU(t)]= g{t,r)}(r)dr, 



(31) 



where g(t, r) is the so-called Green's function in the 
generalized sense. 

To prove this, we consider functions g(t, r) having 
the following properties : 

(P 1) g(t, t) is a continuous function of t, r for 



(P 2) For every arbitrary but fixed t = t {) , a g t ^ 6, 
g(t,T ) possesses n continuous derivatives for 
a St ^ t and r ^t^b such that 



frfrO-O+OjTo) 5^/(7-0-0, Tq) 



dr d/ 

d w -^(r o +0, r () ) d»" 



(»=0,1,. 

r/(r — 0, r () ) 1 



.,w-2) 



(P 3) 8 

aSr ^b, 



dt"- 1 dt 71 ' 1 a (r ( y 

For every arbitrary but fixed t = t in 

k* 

L t [g(t, t )]=-i;^(0^(to) 



for a^t^b but ^r ; 

(P 4) 8 For every arbitrary but fixed t=t in a<r <&, 
Pt(g(t,p ))€Q; 



(P 5) For all functions u M (t) (|*=1, 
the basis (28) 



, k)of 



( fif(«,T) ^(OdtsO 

»y a 



identically in a^r^b. 

It shall be proved that there exists exactly 
one function g(t,r) having these properties which 
satisfies (31). 

Theorem 4: I j there exists a function g{t,r) possess- 
ing properties (P 1) through (P 5), it is uniquely 
determined. 

Proof: Let gi(t, r) and g 2 (t, r) be two functions 
which both possess properties (P 1) through (P 5). 
For an arbitrary but fixed t = t in a<C r< b we define 
the function 

go(t)=gi(t, To) — g 2 (t, Tq). 

Then, according to (P 3) it is seen that 

L[g (t)]=0 for a^t^b, t^r . (32) 

In view of (P 2), g (t) £ C n ~ l ; hence it follows from (32) 
that the jump in the (n-l)-th derivative of g (t) at 
t=r must be zero. Thus, g (t) £ C n ~ l and in view of 
(32), g (t)£K. From (P 4) it follows that g (t)£A; 
therefore, g (t)£KDA. On the other hand, it can 
be concluded from (P 5) that g Q (t) £ (Kt)A) - LL . This is 
possible only if g (t) =0 for a ^t^b. Because r was 
arbitrarily chosen in the open interval a^T^Cb, this 
is equivalent to the statement g x (t, r)=g 2 (t, r) for 
a^t^b, a<C r< b, and in view of (P 1) for the whole 
a^t, r^b. 



8 The subscripts t in Lt[g] and pt(g) denote that the differential operator L and 
the mapping p operate on g(t, t) regarded as a function of t. 
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To state the existence of g(t, r) we first prove the 
following : 

Lemma 1: There exists at least one Junction h(t, r) 
which possesses properties (PI), (P 2), (PS), and 

(PS). 

Proof: For fixed arbitrary r in a<Cr<6, consider 
the ordinary differential equation 



L t [g(t,T)]- 
and the two sets of initial values 
=0 (,=0, 1, 



I] V,(t)v,(r) 



W, r) 
dt v 



n-l) 



and 



b n - l g(t,r) 
dt n ~ l 



<h(r) 



(33) 



(34a) 



(34b) 



Both initial value problems (33)/(34a) and (33)/(34b) 
possess unique solutions hi(t, r) and h 2 (t, r), respec- 
tively. These solutions hi and h 2 are continuous 
functions of t } r in the whole square, a^t,r^b. 
Furthermore, they possess n continuous derivatives 
with respect to t in this square. Accordingly, the 
function 



As(*,T)=- 



{k(tj) 

h 2 (t,r) 



for 
for 



a^t^r 
r^tSb 



obviously possesses properties (P 1), (P 2), and 
(P 3). 

We define the coefficient-functions 



C fi (r) = (h s (t, r), n^t)) 



(m=i, . • . ,£), 



which in view of the continuity properties of h z (t, r) 
are continuous functions of r. Hence, the function 

h(t,T)=h 3 (t,r)-J2c,(r)u,(t) 

evidently satisfies properties (P 1), (P 2), and (P 3). 
Furthermore (P 5) obtains, since 



(hit, t), u ll (t)) = (h B (t, r), u^t)) — c M (r) = 



(M= 



.,*), 



owing to the orthonormality of the u^t). 

Lemma 2: There exists a linear transformation 
which maps every vector x£$l-\-p(K) into one vector 
y £ p(K) such that 

(p- 1 (y),%(0)=o o*=i, . . .,*). (35) 

Proof: Let x£ti-{-p(K) be given. Then there 
exist at least two vectors x x ,x 2 for which 

x=X! + x 2 Xj£Q S™M(zp(K)> 



x l and x 2 are determined to within a vector belonging 
to the intersection illi p(K). With coefficients 



c,= (p~ 1 (x 2 ), u^t)) (m=1, 



.,*), 



the following new vectors are defined : 
k n 

yi=x 1 +Sc M p(^), y 2 =x 2 -2> M p«). 

M=l M=l 

Since 

Sc M p(u M )^Qnp(K), 



then 



x=yi+y 2 y^Q y 2 £p(i£), 



(36) 



and, furthermore, in view of the orthonormality of 
the Up, 

k 

(p _1 (y 2 ) , u,)=( P - 1 (x 2 ) , wj-Sc„(p"VW , if„) 



:(p- 1 (X 2 ), ^)-C M = 0*=1,.., fc). 



(37) 



The two vectors y t and y 2 are uniquely determined 
by conditions (36) and (37). This is rather obvious; 
for any other decomposition x=yi-fy 2 of x which 
satisfies (36) and (37) it follows that p" 1 (y 2 — y 2 )£A[)K. 
However, the relations 



(p" 1 (y 2 -y 2 )^M)=o o*=i, 



.,*) 



obtain only if y 2 =y 2 which then implies yi=yj. 

Hence, it has been shown that we can associate 
with every x£ft+p(l£) exactly one vector y 2 £p(K) 
which fulfills (36) and (37). This was the statement 
of the lemma. 

It is necessary at this point to obtain an extension 
of Green's formula for functions having discontinuity 
properties similar to (P 2). Let the function y(t) 
possess the following properties : 

a ) 2/(0 possesses n continuous derivatives for 
a^t^r, r^t^b, where r is a certain fixed point. 

b) y«( r +0)-y w (r-0)=0 (f=0,1 > . . .,71— 2), 



y K 



) (r+0)- 2 / Cw - 1) (r-0)=c, 



where c is a certain fixed constant. 

For any arbitrary chosen v(t) £ C n , it can be 
shown by repeated partial integration that 

(L[y],v)-(y,L*[v]) = <p(p(y), P (v))-c a (r)v(r). (38) 

This is the generalized Green's formula. 

Theorem 5: There exists a function g(t,r) which 
possesses the five properties (P 1) through (P 5). 

Proof: Let h(t,r) be the function which has prop- 
erties (P 1), (P 2), P 3), and (P 5). For every 
arbitrary but fixed r in a^r^b, insert h(t,r) and 
the basic functions v v (t) (v=\ } . .;>.,&) in the 
generalized Green's formula: 



(L t [h\, fv) = (A, L*[v,])+<p(p t (h),A^))-Mr) 

I Si - 
= <p{p t {h), p(v,))—v,(t) (t=1, 



., **). 
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On the other hand, for a^t^b, 
(P3), 



t^r, according to 



(L t [h], v 9t 



^Mt),V v (t))v P (r)=-v v (T), 



and in view of the orthonormality of the v v (t) 

*(p,(A),p(t>,))=0 (^=1, ... ,**). 
This is equivalent to 

p r (A)6(fi J -Up(iO-L)J-=n+p(Z). 

We now consider the (continuous) linear transforma- 
tion from Sl-\-p(K) into p(K) which exists by lemma 
2. Let h£p(K) with 



( P -(h), «„) = (m=1, 



*) 



(39) 



be the image vector of p t (h) under this transforma- 
tion. Then h is a continuous function of r for 
a^r^b. This is obvious because p t (h) depends 
continuously on r, and the correspondence between 
p t (h) and h is continuous. Let x u . . ., x„bea basis 
of p(K) and 

n 
M = l 

It is clear that the coefficients d M (r) are continuous 
functions of r for a^r^b. Hence, the (uniquely 
determined) function 

h(t, r) = p-\h)=±cUr), P- 1 (X,)^K 

M = l 

is continuous in both variables t and r for a^t,rS b. 
We define the function 

g(t,T)=h(t,T)-h(t,T). 

Then #(£,r) obviously possesses properties (P 1), 
(P 2) , and (P 3) . (P 4) follows from 

p,(0)=Pi(fc)-h£O. 

Finally, as a consequence of (39), 

for fi=l, . . ., k, since A(/,r) fulfills (P 5). 

Thus g(t,r) possesses properties (Pi) through 
(P 5) . This proves the theorem. 

Now we are in a position to prove the following 
final theorem: 

Theorem 6: Let g(t,r) be the uniquely determined 
Junction with the properties (Pi) through (P5) ; 
then 

U(t) = 6[f[t)]=£g(t,r)f(r)dr, 



for every given /(t) £ C°, constitutes the (uniquely de- 
tennined) solution oj (t n ). 

Proof: Let v(f)£A* be an arbitrary function. 
For any fixed r in a<^r<Cb we apply the generalized 
Green's formula (38) to g(t,r) and v(t). With regard 
to (P 4) and p(v) £ 12 - 1 it is seen tliat 

{L t [g],v) = (g, L*[v])-v(r)i 
thus, in view of (P3), 

(L,[g], v)= -Jb,Mv»{T) = {g,L*[v])-v(T) 

for ai^ttkb, a<V<6, t=T. We multiply this 
equation by a continuous function /(r) and integrate 
over r. Because g is a continuous function of t and 
r, it follows that 



-T,M(v»,f)=(G[fh L*[v])-(f,v). (40) 

For/(£) we consider the boundary value problem (29) 
which is solvable no matter how f(t) was chosen. 
Let Uiit) be a solution. For the given function 
v(t) £ A* we obtain 

(L[h 1 ],»)=(/ j »)-S(/a)(^), 



and thus together with (40) 

(£[«,], »)= (fll/],X*W). 

Applying Green's formula (11) to the left-hand 
side we conclude that 

[L*[v], u l ) = (G\f], L*M)-*(p(«i), p(«0), 

and because <p(p(ui), p {v))=Q, then 

((?m-^,i*b])=o. 

This relation is valid for every function v(t) £A*. 
Hence, it follows from the corollary to the alterna- 
tive theorem that 0[J]—u 1 ^C n ; thus, u(f) = G[f} 
£ C n . Furthermore, it follows that 

W]-«J=o, p(0[/]-«O€o, 

which is equivalent to 

L[G[/]]=Z[%]=/-g(/AM, 

p(c7[/])6". 

Finally, it follows from (P5) and (Pi) that 

(«(*), «,.(*))=0 0=1, . . ., ft). 

In applying theorem 6 to the adjoint problem 
(T/f*) we confirm that for &njfi(t) (; C° the uniquely 
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determined solution v(t) of (#n*) is given by 

v(t) = G*[j l ]= P g*(t,r)f(r)dr. (41) 

J a 

We will show that the integral operator G* is the 
adjoint operator of G, i.e.. 



(GU]Ji) = (f,G*\M) 



(42) 



forall/,,/^0 . 

Let v(t) = £*[/,], i.e., 

z*[fl"[/i]]=/i-S(/i,«.)t«. 

pCWiDC^; (fl f *[/J,^=o ( M =i,..,£*). 



Washington, D.C. 



(Paper 64B2-24) 



With another continuous function /(£) we find from 
Green's formula (11), 

(6[f},f 1 )=(6[f],L*[G*[f 1 ]])+j:(f 1 ,u <! )(0lf],n v ) 

v=l 

=(L[G[f]],G*[f 1 ])+<p(p(Glf]),p(G*[f 1 ])) 

+h(fuU„)(G{f},u v ) 



=(/,6 ? *[/i])+Z)(/A)fe^*[/ 1 ]) 

+^(p(g?[/]),p(g*[/,]))+z;(/i,«.) (©[/],«»). 

Since G[fi] and 6r*[/i] are solutions of (it N ) and (*>*), 
respectively, the last three quantities must vanish, 
which leads to (42). 

From (42) follows for the corresponding Green's 
functions, 

g*(t,T)=g( Tf t). 
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